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concentrations from orbit under the assumption of an infinitely 
deep homogeneous chemistry, when the major-element compo- 
sition is simultaneously determined using a gamma ray detector 
[7,8J. Newer calculations reported here show thac interpretation 
of these data can be in error if the hydrogen content of the near- 
surface regolith is not homogeneous. 

Hydrogen is notable for its high efficiency in moderating cosmic- 
ray-induced fast neutrons to thermal energies. The ratio of thermal 
to epithermal neutrons (those in the process of slowing) is a very 
strong function of the hydrogen content. In order to get high 
precision at low H concentration, however, the abundances of 
elements that strongly absorb thermal neutrons such as, e.g., Ti, 
Fe, Sm, and Gd, also need to be determined. 

Recent calculations [7,8] showed that a determination of the 
hydrogen content down to levels of 10 ppm from orbit requires 
only a 1% statistical precision. However, interpretation of these 
data is subject to systematic uncertainties that may be larger than 
the 1% level, depending on the nature and magnitude of H 
abundance variations as a function of depth. Basically, the reduced 
mean-scattering path of thermal neucrons caused by the presence 
of H tends to concentrate neutrons from surrounding material 
into the high-hydrogen-content side of abundance -changing inter- 
faces. For example, comparison of a lunar meteorittc chemistry 
loaded uniformly with 100 ppm of hydrogen down to 200 g/cm 2 
below the surface, and then with 10 ppm H to the 900 g/cm 2 
level, with the same chemistry loaded with 100 or 10 ppm H 
uniformly down to the 900 g/cm 2 level, shows about 2% variation 
in subsurface gamma ray production rates (and hence in thermal 
neutron number density). 

Instrumentation: The neutron detection is made with two 

identical ’He gas proportional counters, one being bare and the 
other being wrapped in Cd. The bare detector is sensitive to both 
thermal and epithermal neutrons, whereas the one wrapped in 
Cd is sensitive only to epithermal neutrons because Cd is a strong 
absorber of thermal neutrons. The abundances of the strong 
neutron absorbers in the soil will be made with a simple scintillator 
gamma ray spectrometer. Although rhis type of detector does not 
offer the high sensitivity of cooled semiconductor detectors, they 
are more than adequate for determination of the abundances of 
Fe and Ti, the principal neutron absorbers on the Moon. Samarium 
and Gd are not easily detected with this approach, but because 
of our extensive database of lunar rock types, these elements can 
be estimated easily from the abundances of other elements that 
will be determined with the GRS, especially K, U, and Th. 

Delivery and Support System: The instruments will be del- 

ivered to the lunar surface via a penetrator with a separable 
afterbody. Using the soil penetration equations of [9] and assuming 
a penetrability comparable to dry silt or clay, a forebody mass of 
30 kg, a diameter of 7.5 cm, and an impact velocity of 100 m/s, 
the forebody will penetrate to a depth of about 10 m while the 
afterbody will remain partially exposed on the surface to maintain 
communication with Earth. Reduction in this maximum depth 
can be readily achieved by reducing the impact speed. Both parts 
will contain neutron and gamma detectors in order to compare 
the differences between the surface and deep abundances. The 
penetrator will be battery operated to provide a life of about one 
week on the lunar surface. 

Conclusions: Calculations of the depth dependence of ther- 

mal, epithermal, and fast neutron fluxes and consequent capture 
gamma ray production rates have shown their utility in deter- 
mining the depth profile of H at the 10 to 100 ppm sensitivity 


level. Orbital surveys of these same data are not capable of such 
a determination. The two techniques therefore complement one 
another. Whereas orbital surveys can provide comprehensive maps 
of suspected hydrogen concentrations due to solar wind implan- 
tation, a penetrator can provide ground truth for a small selection 
of sites to allow estimates of the depth dependence of such deposits 
and consequently an evaluation of the utility of such deposits as 
a source of resources to support human habitation on the Moon. 
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Introduction: Lunar pyroclastic deposits represent one of the 

primary anticipated sources of raw materials for future human 
settlements [1]. These deposits are fine-grained volcanic debris 
layers produced by explosive volcanism contemporaneous with the 
early stages of mare infilling [2,3]. There are several large regional 
pyroclastic units on the Moon (for example, the Aristarchus 
Plateau, Rima Bode, and Sulpicius Gallus formations), and numer- 
ous localized examples, which often occur as dark-halo deposits 
around endogenic craters (such as in the floor of Alphonsus 
Crater). Several regional pyroclastic deposits have been studied 
with spectral reflectance techniques: The Aristarchus Plateau 
materials were found to be a relatively homogeneous blanker of 
iron-rich glasses [4,5]. One such deposit was sampled at the Apollo 
1 7 landing site, and was found to have ferrous oxide and ritanium 
dioxide contents of 12% and 5% respectively [6). While the areal 
extent of these deposits is relatively well defined from orbital 
photographs, their depths have been constrained only by a few 
studies of partially filled impact craters and by imaging radar data 
[7,1]. In this work, we present a model for radar backscatter from 
mantled units applicable to both 70-cm and 12.6-cm wavelength 
radar data. Depth estimates from such radar observations may be 
useful in planning future utilization of lunar pyroclastic deposits. 

Radar Scattering Model: The depth of a pyroclastic layer is 

estimated based upon the ratio of backscattered power between 
a mantled region and an area of unmantled terrain assumed to 
represent the buried substrate. Several conditions are required: 
(I) Only the depolarized (same-sense circular) echo is modeled, 
to avoid consideration of large single-scattering facets; (2) there 
is no volume scattering Within the pyroclastic (i.e. , there are no 
inclusions large with respect to the radar wavelength); and 

(3) the depolarized backscatter originates largely within the re- 
golith substrate, and not at its upper surface. Conditions (2) and 
(3) are supported by previous studies of the eclipse temperatures 
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Fig. 1. 

of pyroclastic deposits and analysis of 70-cm backscattering from 
the lunar regolith [8, 9). 

Input parameters to the model are the radar wavelength, the 
radar incidence angle, and the dielectric parameters of the pyro- 
clastic layer and the regolith substrate. The model used here 
permits analysis of a pyroclastic layer with varying density as a 
function of depth. The electrical properties of an Apollo 17 glass 
sample [6] as a function of density p are used here: 

e, = 2.1 p tan 8 = 0.0037 p 

where e, is the real dielectric constant and tan 8 is the radar loss 
tangent. The real dielectric constant of the lunar regolith has been 
measured to be ~3.0, but may vary from mare to highland regions 
1 10). The incident circularly polarized wave is decomposed into 
horizontal and vertical polarized components, and the Fresnel 
transmission coefficients for each interface are calculated. Atten- 
uation within the pyroclastic is assumed to be exp(-2aL), where 
L is the path length and a is the attenuation coefficient 111) 

a = k(0.5*|(tan5 2 + l) 0 5 - l]) 05 

where k is the radar wavenumber in the pyroclastic deposit. This 
expression is solved numerically in our model for the varying 
density within the lossy medium. 

All the scattered radiation within the regolith layer is assumed 
to be randomly polarized, with some total backscattering efficiency 
S. In our model, we thus combine the H- and V-polarized com- 
ponents within the regolith, multiply by S, and split the remaining 
power evenly between the two polarizations for the return trip 
to the surface. The energy that exits the pyroclastic is assumed 
to have random phase, but may have an elliptical polarization 
due to the differing H and V transmission coefficients. The scatter 
from an unmantled regolith is also calculated, and the two results 
are combined into a ratio of total mantled/unmantled backscatter. 
This cancels the S term in both expressions. The ellipticity of 
the final scattered wave is assumed to be close to unity. The results 
are plotted as a function of power ratio vs. depth for a given set 
of electrical and density parameters. 

In mathematical terms, the H and V backscattered powers from 
a mantled region are 

P h = 0.25 * S *A 2 • (T h 5,p * V ,r + T v !,p ’ T v p,r ) ‘ T h r,p * T h p,s 
P v = 0.25 * S * A 2 ‘ fTY' r * T h p/r + T v 5,p * T/") * T v r,p ’ T v p/s 

where A is the total attenuation factor, and the T terms refer 
to H- and V-polarized Fresnel transmission coefficients for the 
layer interfaces noted in their superscripts; s, p, and r refer to 


space, pyroclastic, and regolith respectively. The backscattered 
powers from an unmantled area are 

Ph = 0.25 ’ S * (T h s,r + T v s/r ) * V s 

Pv = 0.25 * S * (T h s ' r + T/ 1 ) * T v r ' ! 

Results: Figure 1 shows an example of this model for an 

incidence angle of 55°, corresponding to the Earth-based viewing 
geometry for the Aristarchus Plateau. The pyroclastic density was 
assumed to be a uniform 1.5 g/cm ! . The regolith real dielectric 
constant is assumed to be 3.0. Both 70-cm and 12.6-cm radar 
wavelengths are tested, and it is seen that the shorter wavelength 
is attenuated much more rapidly in the pyroclastic, as expected. 

The ratio between average 70-cm returns from the plateau and 
the mean lunar echo is —0.30, implying a depth of 6-7 m based 
on these electrical parameters. This estimate is a lower bound; 
the actual buried highland terrain is likely much brighter than 
the Moon-wide average, implying a greater mantle depth. 

Tests of the above model show that ( 1 ) it is relatively insensitive 
to variations in incidence angle between 20° and 60°, (2) it is 
insensitive to regolith dielectric variations from 2.5 to 4.0, 

(3) the scattered wave ellipticity is >85%, and (4) density vari- 
ations with depth in the pyroclastic, within the narrow range 
permitted by realistic values of 1. 5-2.0 g/cm 3 , are not important 
unless the changes are very rapid. The value of the density at 
the upper surface of the layer is probably a satisfactory estimate 
for the overall density of such shallow deposits. 

Conclusions: This technique may permit remote estimates 

of pyroclastic mantle depths from either Earth-based or lunar- 
orbital radar systems. A 3-km resolution map of the lunar nearside 
is available for 70-cm radar wavelength [12], and ongoing work 
at Arecibo Observatory provides additional data at 1 2.6-cm wave- 
length. We anticipate that such studies will be required in the 
planning stages of a lunar base project for resource assessment. 
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Introduction: When humans return to the Moon (“. . . this 

time to stay ...” [1 1) lunar resources will play an important role 
in the successful deployment and maintenance of the lunar base. 
Previous studies have illustrated the abundance of resource mate- 
rials available on the surface of the Moon, as well as their ready 
accessibility [e.g., 2-5). Particularly worth considering are the lunar 
regional (2000-30,000 km 2 ) pyroclastic deposits scattered about 
the lunar nearside. These 30-50-m-thick deposits are composed 



